Vortex matter in low dimensional systems with proximity induced superconductivity 
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We study theoretically the vortex matter structure in low dimensional (LD) systems with su- 
perconducting order induced by proximity to a bulk superconductor. We analyze the effects of 
microscopic coupling mechanisms between the two systems and the effects of possible mismatch in 
the band structures of these materials on the energy spectrum of vortex-core electrons. The unusual 
structure of vortex cores is discussed in the context of recent tunneling microscopy/spectroscopy 
experiments. 
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I. INTRODUCTION. 



The induced superconducting order attracts consider- 
able interest of both theorists and experimentalists for 
many decades starting from the seminal works on the 
proximity efFectjii^ Recently, one sees a revival of this 
interest in connection with the growing number of exper- 
iments carried out for a variety of new artificial systems 
which include two-dimensional electron gas, graphene, 
semiconducting nanowires and carbon nanotubes, topo- 
logical insulators, etc. Exotic electronic properties of 
these systems''^- can cause quite unusual manifestations 
of the proximity effect. Superconducting characteristics 
of such low-dimensional (LD) systems can differ strongly 
from those in the bulk. Thus the experiments on prox- 
imity induced superconductivity provide a unique possi- 
bility to manipulate the basic properties of the super- 
conducting state. Control of superconducting charac- 
teristics can be realized by changing the doping level 
through the gate potential, which creates, e.g., new types 
of tunable Josephson devices;^ Unconventional gap po- 
tential induces, in turn, unusual quasiparticle (QP) states 
both in homogeneous and in nonuniform superconduct- 
ing phases. For LD systems with a nontrivial topological 
structure one can possibly realize the QP modes with spe- 
cific symmetries of the electron and hole wave functions 
at the Fermi level that describe the so-called Majorana 
fermions in condensed matter i^ii^ 

A standard way of studying the QP states in systems 
with a complicated superconducting order is to look at 
the effects of applied magnetic field on the structure of 
the mixed state. For example, if the bulk electrode is a 
type-II superconductor (SC) one can study the structure 
of vortex lines penetrating the electrode and threading 
also the LD system (see Fig. [Ij. It is the goal of this 
paper to review the basic properties of the vortex matter 
formed in the LD layer. Similar problem of the vortex 
matter in the proximity layers naturally arises when one 
faces the challenge of interpreting the scanning tunnel- 
ing microscopy/spectroscopy (STM/STS) measurements 
in superconductors. Probing the energy and spatial de- 
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FIG. 1: (Color online) Sketch of the 2D layer with multiple 
core vortex structure induced by a bulk type II superconduc- 
tor in the vortex state. Two scales of the induced vortex 
are schematically depicted by orange disks in 2D layer and 
cylinder in the bulk superconductor. 



pendencies of the local density of states (LDOS) by 
STM/STS"'^^ provides information of the spectrum and 
of the wave functions in the superconducting state. An 
important part of this information refers to the structure 
of subgap QP states in the magnetic field bound to the 
vortex core which are known as the Caroli-de Gennes- 
Matricon (CdGM) states^^. A fingerprint of these states 
is the so-called zero-bias anomaly (ZBA)ii seen in the 
STM measurements. Obviously, the intrinsic character- 
istics of the bound core states can be masked or even hid- 
den by the presence of a thin defect layer at the surface of 
the bulk SC. In such thin (possibly non-superconducting) 
surface layer, the superconducting coherence is induced 
by proximity to the bulk SC. The masking effect of the 
defect layer is often difficult to distinguish from more 
exotic explanations based, e.g., on the assumptions of 
the superconducting gap anisotropy (see^^ii^ and refer- 
ences therein) and multi-component structure of the or- 
der parameter-i^ii^. Despite all its simplicity, the model 
assuming the presence of a defect layer at the sample sur- 
face can explain quite a variety of features in the vortex 
LDOS experimental data and provides an instructive ex- 
ample of the vortex matter in the LD systems with the 
induced superconducting order. 

Instead of considering various phenomenological mod- 




FIG. 2: (Color online) Matching of Fermi surfaces in 2D layer 
and in the bulk superconductor in the coherent tunneling case. 
In the simple case of isotropic Fermi surfaces the in-plane 
projections of 3D Fermi momenta Q± coincide with the Fermi 
momentum in 2D layer k2F- 



els of the induced gap potential, in our studies of the 
vortex matter we rather use the general microscopic ap- 
proach developed in Ref. [13 and focus on the physical 
mechanisms responsible for formation of the particular 
gap potential and its symmetry. These mechanisms are 
mostly determined by the nature of the electron trans- 
fer between the two-dimensional (2D) proximity system 
and the bulk SC. This transfer is strongly affected by 
both the mismatch of the band structures in the cou- 
pled subsystems and by disorder in the barrier between 
them. Without disorder and neglecting the band struc- 
ture effects one arrives at the coherent tunneling model 
according to which the in-plane projection of the elec- 
tron momentum is conserved in course of tunneling. The 
induced gap potential is determined by matching of the 
2D Fermi surface with the in-plane projection of the 3D 
Fermi surface (see Fig. [2]) . A generalization of the above 
model can include umklapp processes accounting for the 
Bloch - type single-electron wave functions in both sub- 
systems. In the latter case, the momentum of tunneling 
electrons is conserved only up to certain vectors of the re- 
ciprocal lattices. One more limiting case is the so-called 
incoherent tunneling model which assumes a strong dis- 
order in the tunneling barrier and allows for an arbitrary 
random change in the momenta of tunneling electrons. 
The systematic analysis of these three tunneling models 
shows that the gap potential strongly depends on the de- 
gree of disorder as well as on the band structure effects. 
Based on these models we consider several fundamen- 
tal properties of the vortex matter in the systems with 
induced superconducting order. First, the proximity in- 
duced superconducting gap A2D is responsible for ap- 
pearance of a new length scale in the vortex structure, 
the 2D coherence length, ^2D = Hv^f/^^d or ^2D — 
yJhD2Dl ^2D for clean or dirty limits, respectively. Here 
viF and -D2D are the Fermi velocity and diffusion con- 
stant in the 2D layer. The energy gap A2D depends 



on the tunneling rate F— "— ; for example, A2_d ~ T 
for r ^ A. Since A2D ^ A the coherence length ^2D 
usually is much longer than the coherence le ngth in the 
bulk SC, is = WfII^ for clean or ^s = y/hDg/A for 
dirty limit, where A, Vf and Ds are the gap, the Fermi 
velocity and diffusion constant in the superconducting 
electrode. As a result, all the effects associated with 
overlapping of neighboring vortex cores as well as the 
normal QP scattering at the boundary of the 2D sys- 
tem become much more pronounced than in the primary 
superconducting electrode. There appears, e.g., an in- 
triguing possibility to get a new type of vortex matter 
strongly bonded by the intervortex QP tunneling even 
for magnetic fields well below the upper critical field of 
the bulk superconductor. 

Second, hybridization of the localized QP states inside 
much larger induced vortex cores with the core states of 
primary vortices in the bulk electrode leads to peculiar 
structure of the subgap energy branches. For coherent 
tunneling, the electronic spectrum of a singly quantized 
vortex consists of two anomalous branches crossing zero 
of energy as functions of the impact parameter b. One 
branch, ei(6), qualitatively follows the usual CdGM spec- 
trum eo(6) of the primary vortex; it extends above the 
induced gap where it turns into a scattering resonance. 
The other branch, £2(^)1 lies below the induced gap and 
resembles the CdGM spectrum for a vortex with a much 
larger core radius ~ £,20- Thus, the proximity induced 
vortex in a ballistic 2D layer has a "multiple core" struc- 
ture characterized by the two length scales, 1^5 and £,20- 
Such a two-scale feature does not appear if the proximity 
vortex states are induced by a primary vortex pinned at 
a large-size hollow cylinder tq > ^s, see Refsi2ii2^. 

The spatial and energy dependence of the LDOS in- 
side the multiple core reveals a rich behavior which de- 
pends on many parameters and on the degree of disorder 
both inside the bulk electrode and inside the 2D layer, 
as well as by the barrier disorder. The barrier disorder 
suppresses the influence of the primary CdGM spectral 
branch and leads to broadening of the lower anomalous 
branch e2{b) due to the momentum uncertainty. Impu- 
rity scattering in the bulk and/or inside the 2D layer 
causes further smearing of the spectral characteristics of 
the core states which then approach the usual dirty-SC 
LDOS scaled with the corresponding coherence lengths 
£20- 

And finally, both the nontrivial topological properties 
of the normal state wave functions and the induced pair- 
ing symmetry can affect the presence of the zero energy 
states in the QP spectrum of vortices. This phenomenon 
arises from the wave function symmetry under preces- 
sion of the subgap QP trajectories inside the vortex core 
through the corresponding change in the Bohr - Som- 
merfeld quantization rule for the angular momentum. 

The paper is organized as follows. In section [II] we 
introduce the basic model used further for the analysis 
of the induced superconductivity. The derivation of self 
energies of 2D quasiclassical Eilenberger equations in a 
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FIG. 3; (Color online) 2D normal metallic layer {Z — 0) 
coupled to a bulk superconducting half-space Z > through a 
tunneling barrier. The electron waves depicted by red arrows 
tunnel from the source placed in 2D layer (red solid circle). If 
the energy is smaller than the superconducting gap, they do 
not penetrate deep into the bulk superconductor but undergo 
Andreev reflection to the hole waves (blue arrows) and return 
to the 2D layer. 



vortex state of the bulk superconductor is given in sec- 
tion Hill In section HV] we discuss the method used for the 
calculation of the subgap state structure in the induced 
vortex core. The main results are presented in sections IVl 
and I VII In particular, section |V] contains the results for 
the subgap spectrum and the local density of states in a 
induced vortex state of 2D layer. In section IVIII we dis- 
cuss implications of our analysis for induced vortex core 
states in graphene. We also discuss some further impli- 
cations of a large value of the induced coherence length 
S,2D for the spectral and spatial characteristics of various 
vortex configurations. Some details of our calculations 
are given in Appendix. 



while r is a two-dimensional vector in the normal layer, 
respectively; r is an imaginary time variable in the stan- 
dard Matsubara technique. The chemical potential Ep 
is supposed to be equal in the subsystems. The single- 
particle Hamiltonian in the 2D layer €20 includes the 
kinetic energy and, in general, the lattice potential cor- 
responding to the crystal structure of the normal system. 
For simplicity we neglect the band structure of the bulk 
superconductor. This approximation should be valid for 
a wide class of heterostructures where the Fermi surface 
in the bulk SC is large compared with that in the 2D 
layer. We assume that tunneling is spin-independent and 
occurs locally in time and in space, i.e., from the point 
near the interface R = {r, Z = 0) on the superconduc- 
tor side into the point r in the layer and back with the 
amplitude t{r) that depends on the coordinate of the tun- 
neling center on the interface. Since the tunneling am- 
plitude accounts for certain region of an atomic size in 
the vicinity of tunneling center, the wave function mag- 
nitude at Z = should be considered as an average value 
near the exact boundary of the superconducting region. 
The tunneling amplitude is assumed small in the atomic 
scale. More detailed restrictions for the value of tunnel- 
ing amplitude will be considered later. The tunneling 
Hamiltonian has the form 

+r(r)a+(x)*,(x)] d\ (3) 

where the wave functions in the superconductor are taken 
at the space-time point x at the interface Z = 0. 
The Matsubara Green functions take the form: 



II. 



MODEL. 



Consider a 2D normal metallic layer (Z = 0) placed 
in a tunneling contact with a bulk superconducting half- 
space Z > with a thin insulating barrier between them, 
see Fig. |31 The Hamiltonian of our system has the form 
H = Hs + H2D + Ht, where 



Hs 



d^R 



^#+(X)(e3D-i?F)*.(X) + 



A(R)*+(X)^+(X) + A*(R)*;(X)^I't(X) (1) 



is the part describing the superconductor with the s-wave 
order parameter A(R), csd is the kinetic energy opera- 
tor, and 

H2D = d fd\J2dtix)[e2D-EF]a,i^) (2) 



is the 2D layer Hamiltonian. We introduce space-time 
variables X = (R, r) and x = (r, t) where R is a three- 
dimensional vector in the bulk superconducting region 



and 



(r^aQ(xi)a^(x2)) = (5a/^G(xi,X2) 

(r,#„(Xi)a+(x2)) =5o/jGt(Xi,X2) 

(r,*„(Xi)*+(X2)) = J„0Gs(Xi,X2) 



{Traa{xi)ap{x2)) = iCT^^JF(xi,X2) 

(T,*„(Xi)a^(x2)) = zai^ji^T(Xi,X2) 
(T,*„(Xi)*^(X2)) - iaSF5(Xi,X2) 



(4a) 
(4b) 
(4c) 

(5a) 
(5b) 
(5c) 



etc. Equations for the Green functions can be more con- 
veniently written in the frequency representation a;„ = 
{2ti + 1)ttT. We denote t = ri — r2 and write 

G(ri,r2)= / G(ri,r2;T)e'""^/'"'dr , 
Jo 

skipping for simplicity the subscript. We introduce also 
the Nambu matrices for Hamiltonian and Green func- 
tions 



ho - Ef -A(R) 



^^-l A*(R) esD-EpJ^ 



G = 



G F 

-F+ G 



and denote the inverse operators 



A. Tunneling with umklapp processes. 



G'2d(i") = -^T's^n + To (8) [e2D ~ Ep] , 

in the superconductor and 2D layer, respectively. 

Equations for the mixed Green functions Gt (Ri,r2) 
can be written in the form 

G^HROGtIRi, r2) + dt(Rii)G(Rii, r2)(5(Zi) - 



where Zi > 0, Ri = (Ri±, Zi) and 

t{r) 



t(r) 



t*{r) 



Neglecting the back-action of a thin 2D layer on the 
superconductor, we assume that the superconducting 
Green function Gs(Ri, R2) is a non- interacting function 
that satisfies 

G5i(Ri)Gs(Ri,R2) = iM(Ri-R2) (6) 

in the range Zi_2 > 0. The boundary conditions for 
Gs (Ri,R2) a,t Z — depend on the particular interface 
in the absence of tunneling. This gives 



GT(Ri,r2) = --y G5(Ri,r')t(r')G'(r',r2)dV (7) 

Equations for the Green functions in the layer can be 
written as 



G2-^(ri)G(ri,r2)+t*(ri)GT(ri,r2) 

= id-^hS{ri -r2) 

Using Eq. (O we find 

G2-^(ri)G(ri, r2) - J ET(ri, r')G(r', ra) d\' 

= ihd-^6{ri - r2) , (8) 
where 

ST(ri,r') = (^ _^^t ll ) = ^i*{r,)GUriyW). 

(9) 
One can introduce the momentum representation of 
the Green function^^ 

Gs(Rr,R2) =/00Gs(Qr,Q.)e^'^^--^^^-^ . 

(10) 
and the tunneling coefficients: £(r) — J ,^ ?j f(q)e^'^'". 
The Fourier representation for the Green functions in the 
2D layer is 



The crystal structure of the 2D layer accounts for an 
atomic-scale periodic potential in Eq. ([8]) which mixes 
the Fourier harmonics with the momenta shifted by the 
reciprocal lattice vectors b. Using the Bloch functions 

b 

diagonalizing the single-particle energy operator inside 
the layer 

e2Z3(r)'0m(k, r) = e„^{k)^p„^{k,r) 

one can conveniently introduce the field operators aa,m,k 

aa(i") =y2 77r^aa,,„,kV'm(k,r) . 
^J (2^)' 

The index m enumerates the energy bands. 
Introducing the corresponding Green functions 



{TT0.a,miMAp,m2M'^ — ^aP^mi,m2{'^l7^2) , (12a) 

y) 



{Traa,mi,iii'^P,m2,~ii2) = ^<^aB Pmi,m20^1 7^2) , (12b) 



one can diagonalize the operator G2D i'^ ^'i- ® i^ ^^'^ 
Bloch representation. 



G'2D.m(k) ^-ihfsUn 



e™(k) ~Ef 

e™(-k) - Ei 



(13) 
We assume in what follows that the amplitude A.ind of 
the induced superconducting gap A2D is small compared 
to the interband distance e™ — em' and neglect the in- 
terband scattering. Hereafter we omit the subscripts m. 
At the same time, the transformation from the momen- 
tum to the quasimomentum representation results in the 
mixing of Fourier harmonics in the self energy in Eq. 
([H). Finally, Eq. ^ for the Green functions ^^ takes 
the form: 

G-^(ki)G(ki, ks) - J ET(ki, k')G(k', k2)d^k' 

= M5(ki - k2) , (14) 
with 

ST(ki,k') = ^ /"i+(ki,Qi)G^(Q,Q')4(Ql,k')rf=^Qd'Q' 



4(Q, k) = ^ Wk+bi'(Q± - k - b) 



(15) 



and i^(k,Q^) = 4*(Qj.,k). Here Q = (Q±,g.). The 
above expression for the tunneling coefficients tf, de- 
scribes in fact the umklapp processes caused by the pe- 
riodic crystal potential in the 2D layer. 



B. Coherent tunneling 



in the equation for the Green function becomes 



The simplest model of tunneling assumes that the in- 
plane momentum projection of electrons is conserved 
during the tunneling process: i(Q± — k) = {■5{Q± — k). 
This is equivalent to the assumption that the tunneling 
amplitude t(r) is independent of the coordinate along the 
SC/2D interface. Of course, the quasimomentum con- 
servation is not exact in the presence of energy bands 
since the tunneling mixes the quasimomentum values 
which differ by a reciprocal lattice vector: ^(Qik) = 
iJ2h "k+b5(Q± — k — b). Neglecting umklapp processes 
for simplicity we find from Eq. (J14p 



dQ, dQ', 

(27r)2 



dt^ r 

(ki,k') = — yG°(ki,Q.;k',g'J 



From now on we will use the quasiclassical approxima- 
tion for the Green functions. In order to derive the Eilen- 
berger equations in the 2D layer we follow the standard 
procedure described, e.g., in Ref. [2J. First of all we in- 
troduce the average k = (ki+k2)/2, Q^ = {Qiz + Q2z)/'2. 
and relative k_ = ki — k2, Qz = Qiz — Q2z momenta and 
denote G(ki,k2) = ^(k, k_), Gs{^i,Qiz;'k2,Q_2z) = 
QsCkjQz^k-^qz). Next we apply the operator Gj^j to 
the Green function t/(k, k_) from the right and subtract 
this equation from Eq. p^ . We now transform to the 
quasiclassical Green functions by integrating the result- 
ing equation over d^2 where ^2 = e2D(k)— -Bf- The Green 
functions are to be taken in the vicinity of the Fermi 
surface. Therefore, in the mixed momentum-coordinate 
representation. 



as(k,Q,;r,Z) 
we can put 



^(k,k_)e^- 



gs(k,Qz;k-,qz)e 



k-r+ig^Z Q K-dQz 



(27r)3 



gs{KQz;r,Z) = 5s(k,Q3;r,Z)7ri(5A(6) , 
^(k,r) = g(k,r)^*5A(6) 

Here the standard quasiclassical Green functions are 



g(k2F,r) - - / d6^(k,r), 



5s(Kf,R) 



1 



di^CsiQ,^) 



(16) 
(17) 



^3 = £s(Q) ~ Ep is the normal QP spectrum in the 3D 
half-space, and ^a(C2,3) is a delta function broadened at 
the gap energy scale A. 

At the next step of derivation we note that, in the 
mixed representation, the term 



^/sT(ki,k')G(k',k2)^ 

in J \^T^) 



nidt^ f dQ 



2tt 



gsiQ; r, 0)5(k2F, r) x (5A[e3i)(Q)-^F] 



where Q — {ii.2F,Qz) has the in-plane projection coin- 
ciding with the 2D Fermi momentum k2F- 

Finally, we obtain the quasiclassical Eilenberger equa- 
tion for retarded (advanced) Green functions 

-iftv2FVg(k2F,r) - e[f3g(k2F,r) -g(k2F,r)T3] 

- [ST5(k2F, r) - .g(k2F, OSt] = , (18) 

where ?iv2F = <9e2F)(k)/9k is the 2D layer Fermi velocity. 
For isotropic Fermi surfaces in both the supercon- 
ductor e3_D(Q) = h^Q"^ /2m and the 2D layer em(k) = 
h'^k'^/2m2D, the self energy takes the form 

ST(k2F,r) = yL9s(Q+;r,0)+gs(Q-;r,0)], (19) 
with the tunnehng rate 

V^de 5a [es(k2F, Q.) - Ef] dQ, . 
Jo 

The 3D momentum Q± = (k2F, ^Qsz) lies on the Fermi 
surface of the bulk SC, k'^p+Ql^ = Kp- Provided the 2D 
Fermi surface is smaller than the extremal cross section 
of the 3D Fermi surface, i.e., A;2F < Kp the expression for 
the tunneling rate reads: F = dmP/Q^z- For large 2D 
Fermi surfaces fc2F > Kp the self- energy term vanishes, 
and the coherent tunneling is impossible. The case of 
momenta A:2F — Kp deserves special consideration which 
should take account of a finite delta function width: F ~ 
di2(m/A)i/2. 

The umklapp processes should, of course, modify the 
self - energy part resulting in additional contributions: 



2^(0)/ 



ST(k2F, r) = _^ |uk,^+brsV"'(k2F + b, r) 



where S^^(k2F, r) is given by the Eq. (fT9| 



(20) 



C. Incoherent tunneling 

The coherent tunneling model in many cases oversim- 
plifies the realistic experimental situation. The momen- 
tum conservation is violated, for example, by the pres- 
ence of disorder at the interface. Here we consider an op- 
posite limit of strong disorder, which is sometimes called 
the incoherent tunneling model. This model assumes a 
random tunneling process of electrons through the barrier 
in a way similar to the standard theory of dirty metals 



within the Born approximation-^. We assume that the 
ensemble average of tunnehng ampUtudes is 



i(ri)i(r2) - t's„,5(ri - rs) , 



(21) 



where Sa is the correlated area of the order of atomic 
scale. Following the standard diagrammatic procedure 
we expand the solution for the ensemble averaged Green 
function in a series in the scattering field and split the 
multiple correlators of the t{r) values in a product of the 
above pair correlators. Finally, after averaging the self 
energy ^ becomes: 

ST(ri,r2) ^t^dsaGs{ri,ri;0)S{ri - ra) 

= t^dsaini^siO) (5s(Q; r, 0)) <5(ri - rs) . (22) 

Here 1^3(0) is the normal density of states in the bulk 
material. Angular brackets denote averaging over three- 
dimensional momentum directions. Within the quasiclas- 
sical approach, the resulting self energy to be used in the 
Eilenberger equation (J18p is given by 



tT{r)^ir{gs{Q;r,0)) 



(23) 



where the tunneling rate is F = TTi'3{0)dsat^- This ap- 
proximation coincides with that used in Ref. llTl The 
tunneling rate F ~ t'^/Ep can be expressedii in terms 
of the normal-state tunnel conductance G = 1/ RS per 
unit contact area, F = G / {AttGqV2) ^ EpRo/R , with 
the conductance quantum Gq = e^/irh and the normal 
2D density of states (DOS) 1^2 = rn2D/2TTh^. Therefore 
T/Ep ^ 1 if the total tunnel resistance R is much larger 
than the Sharvin resistance Rq — {NGo)~^ for an ideal 
TV-mode contact with the contact area S. Nevertheless, 
there is a room for the condition F ^ A to be fulfilled 
even for the large contact resistance R^ Rq. 



D. Adiabatic approximation. Range of validity. 

The above microscopic analysis allows us to comment 
on the simplest phenomenological model which is often 
used for description of the proximity induced supercon- 
ductivity, see for examplc )^^'^^" — . Within this model, the 
Bogoliubov ~ de Gennes equations inside the proximity 
superconductor include a phenomenological gap function 
which is postulated to be proportional to the gap func- 
tion A inside the superconducting electrode. Our ap- 
proach shows that this is generally not the case. The 
true equation (jlSp includes self-energies which are com- 
plicated functions of energy, coordinates, and momen- 
tum. In fact, the effective gap function resembles that 
in a usual superconductor only if the bulk SC is homo- 
geneous in space. In this case, the quasiclassical Green 
function is 



5, 



RiA) ^ ^ 



1 



V^^^^W 



e A 
^A* -e 



In this case the self energy is Et = i^gs, for both co- 
herent and incoherent tunneling models. This expression 
also holds if the superconducting gap is a slow function of 
coordinates on distances of the order of ^5. For |e| < |A| 
the self-energy has the form 



ET(r) = 



VWrW 



e A(r 
-A*(r) -e 



(24) 



Only for low-transparency tunnel contact, F ^ A, this 
self energy is nearly off-diagonal on the scale e ^ F and 
can be regarded as an energy-independent effective gap 
function 



iFfae^ 



(25) 



where (j) is the phase of the superconducting order pa- 
rameter. Note that the resulting induced gap does not 
at all depend on the gap magnitude |A| in the bulk. If 
the transparency is finite, the electronic spectrum in the 
induced superconductor has a gap A2D which is deter- 
mined by the conditioniiii 



(e + Si)2-S2=o, e = A2D 



(26) 



Of course, the adiabatic approximation also breaks down 
if the order parameter A varies as a function of coordi- 
nates at distances of the order of coherence length in the 
superconducting electrode, when the self-energies are no 
longer determined by Eq. pi|. 



III. VORTEX POTENTIALS AND GREEN 
FUNCTIONS FOR CLEAN SYSTEMS 

The quasiclassical Green functions in the 2D layer sat- 
isfy the Eilenberger equations (fT8|) . In components, 

-zriV2FV/-2(e + Ei)/ + 2I]25 = 0, (27a) 

ifiv2FV/t - 2 (e + Ei)/t + 22^5 = 0, (27b) 

-ihv2FVg + S2/t - st/ = 0. (27c) 

and the normalization condition g^ — ff^ = 1 with the 
self energies Eqs. ([T^ or ([25]) as effective potentials. 

In this and the following Section we consider the case 
of isotropic Fermi surfaces. Modifications due to the 
anisotropy of the spectrum are discussed in Section IV Bl 
QPs in clean systems are conveniently described by the 
coordinates along their trajectories (see Fig.|4]). A quasi- 
classical trajectory is parameterized by its angle a with 
the X axis, the impact parameter b = psin(0 — a) and 
the coordinate s = pcos{(j) — a) along the trajectory. We 
introduce the symmetric and antisymmetric parts of the 
Green functions as it was done in^^ii^: 

/ = - [C(s) + i9{s)] exp(m) (28a) 

f'^[C{s)-i9{s)]eM~ia), (28b) 




FIG. 4: (Color online) The coordinate frame near the multiple 
vortex core. Primary (induced) core is shown by the white 
(gray) circle. The QP trajectory with an impact parameter 6 
(line AB) passes through the point (p, (j>) shown by the black 
dot. 



where ({s) = C{—s), and 9{s) = —9{—s). The normal- 
ization condition requires g"^ + 9^ + ("^ = 1. Eilenberger 
equations ()27|) can be rewritten in the form 



(29a) 
(29b) 
(29c) 



fiV2F-f- + 2 (e + Si) 6* - 2igEfl == 0, 
as 

hv2F-r - 2 (e + Si) C - 2igY.i = 0, 
as 

fiV2F^ + '^iC^R + '^iOY.i = 0, 
as 



for long distances (p ^ Cs) is independent of the partic- 
ular tunneling model and of the disorder in the bulk SC: 
Si « 0, S2 « Te"!', i.e., S7?, w Ts/p, S/ « Tb/p. How- 
ever, the induced vortex potentials close to the primary 
vortex core are very sensitive to the impurity concen- 
tration and momentum exchange during the tunneling 
process. 

For clean bulk SC, the Green function can be parame- 
terized similar to p8)) with f ~^ fs, C -^ Cs, and 9 ^ 9s- 
The Eilenberger equations have the form of Eqs. (P^ 
with V2F — > I'll = VfCOSXp where Xp is the polar angle 
of the momentum, while Si = 0, S2 ^- A = Ao{p)e^'^, 
and S2 — >■ A*. For energies e <^ Aqo and distances s of 
the order or less than the core size, the functions gs and 
fs are given in Refs i^'^'^^ . 



Cs 



hv\ 



~K 



(e-^-^)Csds', (31) 
P 



2A [e - eo ± i5] ' " hv 

eo{b) = 6A-1 / [Ao/p]e-^(^) ds, (32) 

b\ 



A = / e-^(^) ds ; K{s) = — / Ao(p') dp' 



(33) 



For larger distances, s :^ £_s, the function (s assumes 
its asymptotic expression Cs ~b/p corresponding 

to the boundary conditions Eq. (P5|) . 



A. Vortex potentials for coherent tunneling 



where 



2S;,^=(S2e-*" 
2zS/ = ( Ssic" 



S^e*" 



Sje*" 



The vortex potentials induced in the 2D layer cru- 
cially depend on the tunneling mechanism. For example, 
(oUa) within the coherent tunneling model we get 

(^°^) Si = iTgsiQ, r), S2 - iTfs{Q, r) 



In the present paper we consider the limit of low tun- 
neling rate T ^ A which leads to a small induced gap^i 
A2D = r and long coherence length ^20 ^ £,s- We 
consider an isolated vortex line oriented along the Z axis 
perpendicular to the SC/2D interface and choose the gap 
function inside the bulk SC in the form A = Ao(p)e*'^, 
where (p, </>) are the cylindrical coordinates; Ao(/9) ap- 
proaches the bulk value Aqo far from the vortex core. 
The self energies in the 2D layer, Eqs. p^ or (^5]) . have 
parts with sharp peaks localized at small distances p ^ ^s 
and the adiabatic long-distance "vortex potential" tail 
A2D ^ Fe*"^ at p > ^5 according to Eq. (^5)) . 

In the case of clean bulk SC we use the condition of 
specular reflection at the interface. This can be applied 
for both coherent and incoherent tunneling models since 
any possible disorder in tunneling affects only a tiny frac- 
tion of bulk electrons whose vast majority reflects with- 
out tunneling. For specular reflection, one can use the 
bulk quasiclassical Green functions obtained for an infi- 
nite space. For energies e ^ Aqo, the self energy, Eq. (pS)) 



in terms of the infinite-space Green functions, since 
gsi+Qsz) — gsi^Qsz) for specular reflection. For en- 
ergies e ^ Aoo and distances s of the order or less than 
the core size ^5, we find from Eq. 



(34a) 
(34b) 
(34c) 



Si = 


= -rC5, 


S2 = 


= r[9s ~ iCs]e 


St. 


-- V{9s + KsV 



where C,s and Qs are given by Eqs. ([5T |) - ([55)) . 



B. Vortex potentials for incoherent tunneling 

For incoherent tunneling, we find Si = zF (gs), S2 = 
iV (fs) where averaging over the 3D momentum direc- 
tion is equivalent to the ensemble averaging. To calculate 
the angular average one can separate the Green functions 



8 



into the principal- value part and the delta-functional con- 
tribution. For example, 



9^' 



^C,«(■^) 



ihvi 



-K 



2A[ 



eo] 



± 



■nhvne 



-K 



2A 



-<5(e-eo) 



(35) 

Performing averaging over the polar Xp ^-nd azimuthal 
a angles we take into account the symmetry of the func- 
tions under the s-inversion transformation. As a result, 
we obtain 



Si = 


= -r(Cs(s)) , 




(36) 


Sae-^* = 


= Y\e'^ = ^ad + S^°^ , 




(37) 


^ad{p) -- 


-- p(r/(s)sign(s)/2A[e- 


-eo]) . 


(38) 


We put here 









I{s) ^ 2 (\e - Aab/p)e-^'^''^ ds' . 
Jo 

The off-diagonal components of induced potential are 
split into the localized and the long-range parts, Y^if^ and 
"Sad, respectively. The long-range function "Ead can be 
regarded as an adiabatic induced superconducting gap, 
Sad -> r for p ::!> ^5 and T,ad — >■ for p — ^ 0. Averaging 
over the azimuthal trajectory angle a we find: 



^T,'"" 



31]^°'= 




1-3?- 



V^2^7V 



ehvue ^ 



2Anp^n^p^ - £2 



KEi = -sign(e)r ( ^ 



hvne ^ 



2A^e2 _ ^2p2 



hvi\e-^ 



2AVf^V" 



Here the upper (lower) sign corresponds to a retarded 
(advanced) self energy term, 57 = deo/dfe, and we use the 
notation 



1 



{...)smxpdxp 



for the average over the polar angle Xp of the 3D Fermi 
momentum. Note that our calculations are based on 
the first-order approximation in the small parameter b/p. 
According to Eq. ([30l) the symmetrical E/(— s) = S/(s) 
and antisymmetrical Ei^(— s) = — Ei^(s) parts of the off- 
diagonal self energy term 1^2^^^'^ can be rewritten as 
Sfl = Eae-^-^s/p and E/ = Eae^^'^fo/p. 

The self energy obtained above affects the vortex core 
states in 2D layer in two different ways. The adiabatic 
part of the induced vortex potential leads to the Andreev 
localization of QPs with energy smaller than the induced 
superconducting gap F within the induced vortex core 
at distances of the order of £,2D ■ This forms the CdGM 



anomalous branch e2{b) as in an usual superconductor 
with the corresponding maximum intrinsic gap F. An- 
other part of the self energy exponentially decaying at 
p ^ ^s contains information about the CdGM states in 
the bulk SC; it affects the 2D-layer QP behavior at small 
scales. The adiabatic large-scale part of the self energy 
(at jO 3> ^s) is universal; it does not depend on the tun- 
neling models and on possible disorder in the bulk SC, 
while the short-scale induced vortex potential localized 
at small distances does crucially depend on these factors. 
Both terms in the induced self energy form the two-scale 
local DOS (LDOS) radial profile. 



IV. SCALE SEPARATION METHOD 

A natural way to solve Eqs. ([^ is to apply the scale 
separation method. We introduce a distance po satisfying 
^S ^ Po ^ ^2D and consider the Green functions in two 
overlapping spatial intervals: (i) /o < po and (ii) p ^ po- 
Next we match the solutions in different spatial domains. 



A. Large distances 

At low energies e <C Aqo and large distances p ^ £,s 
the induced vortex potential is given by Eq. (|25|) . QPs 
propagating along the trajectories with impact parame- 
ters b > S.S that miss the primary vortex core are affected 
only by this long-distance (^2D 3> ^s) part of the induced 
gap potential. In the low energy limit e < F ^ Aqo the 
appropriate boundary conditions far from the induced 
vortex core (p ^ S,2d) are 



Ts/p 



Vr^ 



C = 



-Tb/p 



VP 



,9 = 



Vr' 



(39) 



For both tunneling models and arbitrary disorder rate 
inside the superconductor and for p^ ^s Eqs. (|29l) take 
the form: 



hv2F^ + 2e6' - 2igTs/p ^ 0, 
ds 

do 

fiV2F-r - 2eC - 2i.gF&/p = 0, 
ds 

dg 



hv2F^ + 2i9Tb/p + 2iCVs/p 
ds 







(40a) 
(40b) 
(40c) 



The functions g and C are even in s while is odd, 
so we can consider only positive s values. We obtain 
the solution of the above equations using the first-order 
perturbation theory in the impact parameter b: w{s) = 
wo{s) + wi{s), where w{s) = {(,9,ig) . As we shall see 
later, this approximation holds for |5| <C C2D- The zero 
order in b solution reads 



wois) 



1 



Vr2 



-uo{s) 



c 



Vr2 - £2 



u_is) , (41) 



where 



M±(s) 



Vr2 - £2 

±e 

±r 



,±As 



, uo(s) = r 



and A = 2\/V^ — e^ /hv2F- This solution satisfies the 
boundary conditions 5 = — ie/Vr^ — e^^ C = and 
6* = r/\/r2 — e2 for s — >• 00 and e^ < T^. The first order 
correction wi can be written as 

^,(,)^^(f)|L + ^M£± + %M£=, (42) 



where 



Vr2 - £2 vr2 - £2 vr2^ 

^2DCo{s)^2Cb e-^'-, 

Js P 

Js P 



^2DC-{s) = -h / e 



,As 



ds 



(43a) 
(43b) 
(43c) 



The lower limit of integration in C_, Scj has to be taken 
as Sc ~ ^s for trajectories that go through the primary 
vortex core, & < CS) so that the logarithmic divergence is 
cut off at the distances ~ ^s where the long-range vortex 
potential Eod ([38]) vanishes. For 6 :^ ^g we have Sc = 0. 
The perturbation approach holds as long as Co ^ C and 
C+ < 1, i.e., as long as |6| < ^20- For s > £,20 the 
coefficient Cq decays faster than exponentially, while 



2Vr2 - £2 p 



so that C approaches —{h/p)T/\/T'^ — e^ and the correc- 
tions to 9 and g vanish as it should be according to ([M]) . 
For a small distance s = so defined as p^ = s^ + b^ we 
have 



C(so)-C + C+(so) + C-(.so), 
1 



0(5o) 



Vr2-£' 



{r-£C-HrCo(so) 



'2 _ g2 

+ £[C+(,so)-C_(so)]}, 



ff(so) = 



Vr^^7 



{-e + rC-£Co(so) 



r[C+(so)-C_(so)]}. 



(44a) 



(44b) 



(44c) 



Expressing the coefficients Co and C-)- in terms of the 
energy 



(h^ 2r26 



hv 



2F 



(45) 



of bound states in the induced vortex core, with rj 
6d/|6|, Co = -2CC+ = C£2(6)/r, we find 



C[£-£2(6)]=r-££2(6)/2r 



(46) 



According to Eq. (1461) £2(&) is the only spectrum branch 
in the energy interval |£| <C Aqo. The Green function is 



ffW 



-le 



iVC 



,-As 



Vr2 



Vr2 



»£C'o(s) 

Vr2Tr^ 



iT 



x/r2 - £2 



[C+(s)e-^''-C_(s)e" 



-Asl 



(47) 



For s » ^213 we have Cq -> 0, C+e^'' - C_e~^'' -)► 0, so 
that the first term is the homogeneous background while 
the rest terms describe the vortex contribution. To ob- 
tain the retarded function for £2 > r2 one has to continue 
\/r2 — £2 analytically throughout the upper half-plane of 
complex £ keeping ^^/T'^ — £2 > 0. 



C. Matching for small impact parameters 

To find the Green functions for small impact parame- 
ters b < £s one has to match Eqs. p4)) with the solution 
obtained in the vortex core region. For small s < sq we 
assume that the even parts of the Green function 17(5) 
and C(s) are nearly constant in the interval < s < sq. 
Integrating Eq. (j29bp over s from to sq along the tra- 
jectory we find the matching condition 



hv' 



'2F . 



{so) 



C(so) 



So rso 

El ds + ig{so) / S/ ds . (48) 
Jo 



Equation pS)) determines the constant C. Its poles de- 
fine the eigenstates of excitations as functions of energy 
and the impact parameter. While deriving the effective 
boundary condition (H51) for b < £s, one needs to sepa- 
rate the exponentially converging parts E'j"^ at s ~ ^s 
from the long-distance, s ^ £,s, asymptotics of 'Sij. For 
£ ^ Aco the long-distance expressions, Eq. (|25|) . yield 
El ^ 0, E/ ^ Tb/p. Therefore, 



B. Matching for large impact parameters 

Far from the primary vortex core at impact parameters 
£s ^ b <^ £,2D the perturbation result Eqs. (HH]) can be 
applied along the entire trajectory so that one can put 
Sq = Sc = 0- The boundary condition for an odd function 
requires 0(0) — 0. Since in this case C_(0) = 0, we find 
from Eq. (|44bl) 

r + £C+(o) = £C-rCo(o) . 



E/ ds — 



' ds 



So 

r / b/pds 

Us 



E^°" ds 



Tbln{so/£s 



(49) 



while /„ " El ds can be extended to infinity. The localized 
self energy parts Ei, Ej°'^ determine the small-distance 
LDOS and spectrum of excitations and depend on the 
particular tunneling mechanism. 
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V. MULTIPLE VORTEX CORE IN THE 

CLEAN LIMIT. QUASIPARTICLE SPECTRUM 

AND DENSITY OF STATES. 

A. Isotropic Fermi surface 

In this section we consider an idealized picture with- 
out any disorder. For large impact parameters, b 3> 
£,s, the corresponding solutions for the Green functions, 
Eq. (|47)) . coincide with the standard CdGM expressions 
where the gap value is replaced with F. The correspond- 
ing anomalous spectrum for 2D excitations is given by 



Eq. 



23,28 



This modified CdGM branch dominates in 



the LDOS at large distances p:^ ^s- 

The normalized FDOS is defined as an average over 
the trajectories: 



7V(r,e) 



2Tr 






^ p^ - 62 TT 



where Nf{s,b) = [g^{s,b) — g^{s,b)] /2, s = pcosa', 
and b = — psina'. For |e| < F, a nonzero LDOS comes 
only from the vortex contribution of the second and third 
terms in (|47|l due to the presence of a pole in the coeffi- 
cient C according to Eq. (|46l) . The Green functions and 
LDOS reach their long-distance values g = —ie/y/T'^ — e^ 
and N = 3fi|e|/Ve2 - F2 as p -> oo. For p > Cs 
the trajectories with large impact parameters b > ^s 
give the main contribution to the LDOS. In the region 
^S -^ P -^ £,2D we get the angle-resolved density of states 
in the form: 



N,is,b) 



VF2-e2(F2-e2/2) 

f2 

X7r(5[e-e2(6)] , 



kl<r (50) 

Vi^^P[F2-e2(6)/2] 

Ne[s,b) = : — ——2- — -— , e > F . 51) 

sign(e)F^[e-e2(6)J 

Thus, the corresponding LDOS in the energy interval 
|e| < F has the only peaks at e = e2{±p): 



N{p,e) 



N.iV7 



db 



'''■"^- 


-62 


VF2-e2(l- 


- e2/2F2) 



= n— — ; ^^ - / "- 2 . (52) 

For energies above the induced gap, |e| > F, for the same 
distances the LDOS is monotonically increasing with |e| 
to its normal state value: 



Nip,e)^Ve 



\e\ , (I~e2/2F2) 

2r2 v^^^^IR 



(53) 



A trajectory with a small impact parameter 6 < fs 
can be divided into the part far from the primary vortex 
core, and the region inside the core. Far from the core 



r 



-r 







,^^^ 


S^)^ 


^Q 


^^(b) 


^^^>^^^^^n 




e2(b) 


b' b* 




^S 


.._—-—■ — 




Imejb) 




1 



FIG. 5: (Color online) Two localized branches, ei(b) and e2(fo) 
of the spectrum, Eq. (|59|) . in the limit of coherent tunneling, 
for e < r. h* is defined as ei(6*) = F — 0, while 6' corresponds 
to5ftei(6') = r+0. The spectrum satisfies £1,2 (—6) = — ei,2(&)- 



the solution is found using the vortex potentials Eq. (|25l) . 
The self energies of the primary vortex in Eq. (jlSp have 
poles at the usual CdGM energy eo(6) with the corre- 
sponding wave functions exponentially localized within 
P ^ Cs and regular parts extending over large distances 
p -> ±oo^i2^: 



Vd. 



-rCs 



(54) 



Note that the localized part Yr^'^ of the effective or- 
der parameter S2 has the coordinate dependence Yi^'^ = 
iY}f'^{b, s)e'" with zero circulation, unlike its adiabatic 
part (HSl) S2(p > "^s) = Te*"^. As we wiU see below it 
is this different angular dependence of the effective gap 
asymptotics, which leads to the formation of a "shadow" 
of the bulk SC anomalous branch in the excitation spec- 
trum and LDOS in the 2D layer. 

Using Eqs. (I43|) for the long-distance part of the tra- 
jectory we find 



Co (so) 



2C6, 1 
7 — ^"^1" 

t,2D Aso 



C+(so)±C_(so)«--^ln-^ 

?2D A^S 



^2(6) 

2F 



(55a) 
(55b) 



We now match the asymptotic solution Eqs. ((31]) ob- 
tained for s > So with the solution for the short-distance 
part of the trajectory, Eqs. (p4|) and (|3T|) - (|33)) . using 
Eq. (gHl) and Eq. (gS])- As a resuh. 



C 



S,2D[e - £2(6)] + 2[F - v/r2 - e2 _ !!i(^] p Q^ ds 

££2(6) 



^2DT + 2ej Corfs-6D-2p 

-(F + Vr^37i)!^^°° Cods' 



where Co(s) is the localized part of C,s and 

hv\\ 
2[£-£o(6)] 



, (56) 



(57) 
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■~\!v[r^p=2.5k 




P=15^c 



-A 



-r r 



FIG. 6: (Color online) LDOS in logarithmic scale for coher- 
ent tunneling in the clean limit. Curves, taken for different 
distances p from the vortex center, are vertically shifted for 
clarity. The peaks in LDOS exist up to distances ~ ^2d ■ Here 
A/r = 5, g„ = 1. 



Here we put g = i(^Q and replace the cutoff parameter in 
dJS]) by 77 = C2d/Cs- For 6 » ^s the contributions from 
the primary vortex core proportional to L Co ds vanish 
since the trajectory misses the core, and Eq. (|56|) goes 
over into Eq. P^ . 

For small b <C S.2D the Green function has a pole when 

Pie,b) ^ [e - e,{b)][e - eoib)] 

+ q, [r2 - rv/r2 - e^ - ee2(6)] = (58) 

where Qv = v\\/v2f- One can show that including the 
higher order terms in the parameter e2{b)/T the corre- 
sponding energy dispersion relation takes the form: 

[e-e2{b)][e-eom 



rqv 



r-v/r2-[e-e2(6)]2=0 (59) 



For ^ ^ Cs ; the cutoff parameter in Eq. P5|) should be 
replaced with 77 = ^2d/£,s- 

The resulting two-scale spectrum is shown in Fig. [5] 
There are two real- valued branches in the range |e| < F 
crossing zero of energy as functions of the impact pa- 
rameter and one complex-valued branch in the range 
F < |e| < Aoo- The lowest-energy branch e2{b) has a 
scale ^2D as a function of the impact parameter: For 
b < ^2_D it is given by Eq. (|45)) with the proper cutoff 
parameter -q as discussed above and saturates at e = F 
for b » ^2D- The branch ei{b) has a scale £,s'- For 
e < F it goes slightly below the CdGM spectrum eo(^) 
of the bulk SC, ei(&) = (1 + qv/2)-ho{b). Above F the 
spectrum transforms into a scattering resonance due to 
the decay into delocalized modes propagating in the 2D 
layer: ei(6) = eo{b) - iTq^ for |e| > F. Since Eq. ((55)) 
determines a pole of the retarded Green function in the 
lower half-plane of complex e, the square root in Eq. ([5^ 
should be analytically continued through the cut going 
from ~oo to — F and from F to +00. As a result, ei(fe) 
has a discontinuity at ei = F with fe'/'Cs ~ 0-29 and 
b*/^s - 0.42. 

The two branches appear due to the presence of two 
sub-systems, the bulk SC and the 2D proximity layer, 



each with its own anomalous branch. The existence 
of two anomalous branches follows also from the index 
theore m^^i'^" . Indeed, its application requires that both 
zero of the quasiclassical Hamiltonian at the Fermi sur- 
face and its singularity at e = eo{b) are taken into account 
when calculating the topological invariant. As a result, 
the number of anomalous branches is increased up to 2 
for a single-quantum vortex. 

The multiple-branch spectrum results in multiple-peak 
structure in the LDOS (Fig.lH]), which appears to be most 
pronounced deeply inside the primary core (at distances 
p < ^s/^2D when ei < F) thus illustrating the two-scale 
structure of the vortex core. The LDOS is obtained from 
the angle-resolved DOS (normalized by its normal state 
value) N^{s, b) = [g^{s, b) — g^{s, b)]/2 averaged over the 
trajectory direction. 

The angle-resolved DOS for small energies |e| <C F and 
p'^£,S reads 

(60) 

Here we neglect the terms ee2{b)/T'^ and e2(fc)/ei(&) and 
put eo(&)/ei(6) — l+gi,/2 according to low energy asymp- 
totics. In this case the LDOS 



N{p,e) = 3? 



Fg^ 



2y^?(^ 



^^k±^ (61) 



reveals a two-peak structure vs energy at e = £1.2 (p)- For 
|e| ~ F, one can neglect £2(6) and obtain: 



[e-eo(&)] F + Vr2 - e2 + g^Fe = . (62) 

For |e| > F the dispersion relation is complex valued and 
for retarded functions takes the form: 



ee- 



eo(&)]+g,,F F + isign(e)\/^ 



F2 







(63) 



The latter equation describes the resonant states in the 
2D vortex core which decay into the QP waves propagat- 
ing in the 2D layer above the induced gap. 

Finally, the whole spectrum structure, shown in Fig. [SI 
has two anomalous branches: (i) one of them 62(6) is 
completely real-valued and follows the CdGM spectrum 
for the superconductor with homogeneous gap F; (ii) an- 
other one is close to the bulk CdGM spectrum, but has 
a discontinuity at e = F, where it becomes essentially 
complex. 

Thus, the LDOS for energies above the induced gap 
|e| > F and smaU distances p,b < S,s reads 



ATf ^ -f^^ ^ 9.r2 



2e 



K 



^g2 _ r2 - iF 



^(e2 + g,F2 + z(7„FVi^^n)2 - e^^l^p) 



(64) 
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2D Fermi surface 



(b) 



FIG. 7: (a) An example of anisotropic Fermi surfaces showing 
a spherical 3D Fermi surface on top of a part of a 2D Fermi 
surface in the layer shifted from the center of its Brillouin 
zone. The closed loops show the 2D Fermi line and its pro- 
jections onto the 3D Fermi surface. The directions of the 3D 
Fermi velocity projection vsf on the plane 2 = in the bulk 
does not coincide with that in the 2D layer, V2f. (b) Differ- 
ent points ri and r2 specified by Si and S2 on a 2D trajectory 
with given impact parameter b belong to trajectories in 3D 
with different impact parameters &i and 62. 



and has the only peak at e = ^ei{p) of the height ^ 
r^/eKp) for p > £,g/^2D- In the opposite limit of rather 
large distances p > (,g/£,2D at |e| > F, the spectrum 
reduces to the CdGM spectrum with a finite broadening: 



ei(&) = eo(&) - «rq„ 



(65) 



The LDOS has a small difference from its normal state 
value iVo = 1: 



7V(p,6) = l + %^3? 



iV 



2e2 ^{e + tq,rr-el{p) 



(66) 



The LDOS in the whole energy range (|61[ IMl) has two 
or even three peaks for such distances. The latter case 
is realized at the distances corresponding to 6' < 6 < 
&*, where the spectrum vs the impact parameter has 3 
anomalous branches. 

The numerical LDOS patterns have been obtained by 
the subsequent solving of the two sets of Eilenberger 
equations in the Riccati parametrization'^^: first, we cal- 
culate the Green functions in the bulk SC using the ap- 
proximation Ao(p) = Aoop/\/p^ + Cs and next we solve 
Eq. (dll) in the 2D layer using Eq. ([T9| . 



one can also apply the method of scale separation in the 
same manner as we did it in Sec. IIVI The consideration 
for the region of large impact parameters does not differ 
significantly, such that the solution for the Green func- 
tions together with the matching conditions look simi- 
lar to Eqs. gl]), gSI), dMI) and gH]), gS]). However, 
the region of small impact parameters of the order of S.s 
gives an essentially different result. The main distinction 
is that the directions of QP trajectories which are de- 
termined by the group velocities de2D/dk and deso/dQ 
for given in-plane momentum in 2D and 3D systems do 
not coincide. Fig. Elja). As a result, the integral in Eq. 
(|48)) along a 2D trajectory involves trajectories with dif- 
ferent impact parameters used to parameterize the 3D 
Green functions, see Fig.ITJb). Within the quasiclassical 
approximation, the integral will thus give an imaginary 
part which comes from the delta function at the 3D core 
spectrum and a real contribution from a smooth depen- 
dence. The spectrum e2{b) at small impact parameters 
thus becomes broadened and shifted from its initial po- 
sition. The imaginary contribution appears due to the 
coupling of the QP trajectory in the 2D layer with a 
quasiclassical continuum of trajectories inside the super- 
conductor corresponding to different impact parameters. 
This coupling results from the non-conservation of the 
angular momentum in the anisotropic system. The imag- 
inary contribution will be also present if the primary-core 
spectrum is broadened by disorder or inelastic scattering. 
The situation is in many respects similar to that for inco- 
herent tunneling model discussed in the following section. 
Of course, attributing the origin of the imaginary part of 
energy for the anisotropic case to the continuum of states 
in the bulk SC, we ignore the angular momentum quan- 
tization in the primary vortex core. The true quantum 
mechanical consideration accounting for the level quanti- 
zation could possibly change this conclusion and lead to 
a real- valued energy spectrum for ideal systems without 
disorder. 

In this Section we consider the low-energy behavior of 
the Green function at small impact parameters, b <^ ^s 
where the spectral energy e2{b) is very small and can 
be neglected. We assume that the trajectories in the 
bulk SC and in the 2D layer do not coincide; the Fermi 
velocities V2f and v^p are at an angle Sa to each other, 
see Fig. [TJ^b) . The impact parameter bs and trajectory 
coordinate ss in the superconductor are coupled with the 
ones in the 2D layer {b and s) through 

bs — P sin((/) ~ a + Sa) — b cos Sa + s sin Sa , 
Ss = P cos{(f> — a + Sa) = s cos Sa — b sin Sa . 



B. Anisotropic Fermi surface 

Here we briefly discuss the effects of anisotropic Fermi 
surfaces in 3D and/or 2D systems. We will be interested 
only in main distinctions which the anisotropy causes 
within the coherent tunneling model as compared to the 
isotropic case considered above. For anisotropic surfaces. 



As we know, at small distances p <S^ S,s the part ^5 of 
anomalous Green function fs in the bulk superconductor 
is large compared with Cs >• ^s- Neglecting the latter, 
we have for the self energies in Eq. ([50]) : 



Efl — Si sinSa, 
Ttj = El cos Sa, 



(67a) 
(67b) 
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The diagonal self energy is Si = iTgs ~ — rCs- Note 
that the self energies depend on trajectory coordinate s 
in 2D layer through the impact parameter bs = b cos Sa+ 
s sin Sa in the bulk superconductor and do not possess 
definite symmetry in s. Therefore, one needs to consider 
the region inside the primary core more carefully allowing 
for contributions from even and odd components of the 
corresponding functions. 

As in Sec. |lV]we use the scale separation method and 
subdivide a 2D layer trajectory with a small impact pa- 
rameter b '^ ^s into the long-distance part far from the 
primary vortex core, and the region inside the core. We 
introduce a distance p' satisfying ^s ^ p' ^ £,2D and 
consider the Green functions in two overlapping spatial 
intervals: (i) p ^ p' and (ii) p ^ p' ■ Next we match the 
solutions in different spatial domains. Far from the core 
the solution is found using the vortex potentials Eq. (1251) . 

In the region inside the primary vortex core the self 
energies play the most important role. Using the ap- 
proximation (j67p . Yjr — Yjisai5a, E/ — YiiCOsSa and 
neglecting e, Eqs. (P^ at small distances s < sq can be 
written in the matrix form 



d 2E1 J 

— -t& -f Aw ~ 

as nv2F 



As in Sec. IIVI we use the vector w 
constant matrix 



(68) 
(C, 0, ig)^. The 






1 


— sin 6a 


-1 





— cos Sa 


sin 5a 


- cos 5a 






A^ 



has threefold degenerated zero eigenvalue, therefore the 
solution of Eq. (j68|) can be written in terms of mutu- 
ally orthogonal eigenvector vq = (— cos (5a, sin (5a, 1) 
and adjoined vectors vi — (— sin5a, — cos5a,0) and 
1)2 = (cos 5a, — sin 5a, 1) : 

w = Cq{s)vo + Ci{s)vi + C2V2 , (69) 

where Avq = 0, Avi = vq, and Av2 = 2t'i. Therefore 



dCijs) ^ , 2Si(g) ^^ 

dx hV2F 

dCo{s) 2Ei(s) 

Ci(s) 



(70a) 
(70b) 



dx hv2F 

The solution is 

Ci = C° + 2C2l{s) , Co = Co" + C?/(5) + C2l{s)\ (71) 
where Si = — F^s and 

2F 



I{s) 



hV2F Jq 



Cs{s')ds' 



(72) 



The three equations (jHSj) arc not independent because 
of the normalization g'^ + (^ + 9'^ = 1. Therefore, the 
three coefficients C", C2, and Cq are coupled through 

(C?)2 = 4C2C0". 



Multiplying Eq. (|69)) by the vectors u^ and v^ and 
using Eq. (|71l) one can obtain 4 conditions at s = zLsq 



(C cos Sa — 6 sin Sa + ig) 



±sa 



C2 



(73a) 



(C sin 5a + 6* cos 5a) ±^^ = -C? - 2C2/(±so) (73b) 

Excluding the coefficients C2 and C" we find 

[(cosSa — 6 sin Sa + ig]^ —0, ('''4a) 

[C sin Sa + 6 cos Sa] ^^ 

+2Iodd {C cos Sa-9 sin Sa + ig}^^ =0, (74b) 

where the integral I{s) — Ieven{s) + Iodd{s) in Eq. ([7^ 
separated into even Ieven{~s) — Ieven{s) and odd 
loddi^s) = -Iodd{s) parts and [x]^^ = x{so) - x{~so) 
and {x} ^ = a;(so) + a;(— so). The integral 



lodd = 



hv 



2F J- 



So 



Cs{s')ds' 



takes the form 

Fi'ii 



lodd — 



-V 



W2F2Asin5a 
exp 



00 — i<'(^ cot (5q — 6 sin (5a) 



2Av2F^ sin 5a 



-K 



[e — eo (6 cos 5a + z)] 
e cot Sa b 



dz 



n 



sin 5a 



(75) 



where we put s sin 5a — z. The second term comes from 
the delta-function contribution at one of the primary 
core states, see Eq. (PT|) : the upper (lower) sign cor- 
responds to retarded (advanced) function. For Sa < e/A 
the second term disappears while the first gives the real 
pole contribution which is equivalent to Eq. (|57p. One 
concludes that the imaginary part disappears only for 
trajectories which are almost parallel (within an angle 
Sa < e/A). For Sa ^ e/A the first (real) term vanishes 
since the expression under the integral becomes odd in 
z. Note that for e = and 6 = the real term vanishes 
exactly. 

Equations (j74p are the matching conditions with the 
solution in the large-distance region, s > sq. They are 
generalizations of the matching condition Eq. p8|) de- 
rived earlier for the isotropic situation. The two condi- 
tions Eqs. (|74p determine the even and odd parts of the 
Green functions. 

The long-distance solution is found in the same way 
as in Sec. IIVI However, it does no longer have a definite 
symmetry with respect to s ^ — s. We separate the even 
and odd components w — Weven + Wodd and consider 
both s > and s < 0. In this Section we only discuss 
the behavior of the Green function for low energies and 
small impact parameter. We thus neglect the corrections 
to w proportional to b/^2D- lu this case Weven is given 
by Eq. (|41l) where now 



Vr2 - e2 

u±{s) — \ ±esign(s) I e 

±F 



±A|s| 



, Uo(s) 
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and 



Wodd 



C'sign(s) 
x/r2 - £2 



u_(s) 



Equation (|74ap gives 



C 



sin Sa{T -Ce) 



\/r2 - e2 cos (5a - r 



(76) 



(77) 



Using Eqs. (PTI) . ([75)) . and ([771) we find the combinations 
C(so) ± C(-so), ^(so) ± d{-so), and 15(50) + igi-so) in 
terms of the coefficient C. Next we insert these combi- 
nations into Eq. (|74bl) and find 



C[e ~ 2ZIodd] = r + 2eIodd cos 5a 



Z = 



[Vr2 - £2 cos 5a - r] 2 - £2 sin2 5a 

\/r2 — £2 — r cos Sa 



(78) 



(79) 



where 



Equations (|78|) , (|79l) are the counterparts of Eq. (l46l) for 
an asymmetric case and transform into it for 5a — >■ 0. 

For £ ^ r we have Z = 1 — cos 6a. For 5a > £/ A the 
integral /odd Eq. ([75]) has only imaginary part. There- 
fore, 



C 



T 



£ — £2(6) ± 17 



where 



_ 7rt>||rtan(5a/2) ^^(^^^ 



Ai;2F^ 



rVA, 



(80) 



(81) 



and po = |fc/ sin 5a I . In Eq. ([SO]) we include the energy 
£2(6) which can be obtained by more detailed calculations 
taking into account corrections due to b/p in the same 
way as in Section ITVl The function K{pq) decays expo- 
nentially as e^P"'^^ for impact parameters larger than 
the primary core size, b> ^s- 

Therefore, the imaginary term in (j80l) does not disap- 
pear unless Sa is very small. It results in smearing of 
the adiabatic energy level e2{b) <C F and in a Lorentzian 
behavior of the DOS due to tunneling into the primary 
vortex core states. We remind that this result is obtained 
within the quasiclassical approximation. 



and, thus, use the results of the previous section for the 
Green functions. The self- energy potentials are now ob- 
tained by averaging the Green functions Eqs. (|3T|) - (|33|) 
over the trajectory direction: E^ — iT{gs). This aver- 
aging does not affect, of course, the induced gap func- 
tion ()25p outside the primary vortex core and, thus, the 
spectrum £2 survives the influence of the tunnel barrier 
disorder at least for b > ^s- On the contrary, the sub- 
gap branches localized within the primary vortex core 
are completely destroyed. Such dramatic consequence of 
the momentum scattering is caused by the averaging of 
electronic wave functions with different impact parame- 
ters and consequent loss of information about the CdGM 
states of the primary vortex. A natural consequence of 
the momentum scattering is the appearance of a finite 
broadening of energy levels for trajectories with small 
impact parameters b < ^s- Matching the solutions in the 
core and at large distances gives the expression for the 
coefScient C for 6 < ^5 and |£| ^ F, 



C 



£2(6) 



2Vr2 



hv- 



'2F 



El ds — 



2T 

hv- 



'2F 



E^°" ds 



2e 

hV2F 



Y}r ds 



(82) 



Since |Ei| - jEf^l - F the pole of the coefficient C is 
located at small energies e < F2/A <C F. Thus, for £ ^ F 
the expression for this coefficient takes the form 



C 



e-£2(6) + 



2D Jo 



(Si 



E^°") 



ds 



(83) 



The localized self energies Ei and Ej°'' can be neglected 
for £ ~ F. They also vanish for |6| ^ ^5. In both 
these limits, Eq. (|S2I) transforms into Eq. (gS]). The 
integral term in the equation above can be written 
in terms of its real /3(6) = (3i{b) — /3i(6) and imaginary 
7(6) — ^i{b) — 71(6) parts as follows: 



2 1'°° 
■— (El - E^-) ds = -m ± Z7(6) 

?2D Jo 



(84) 



Here upper (lower) sign corresponds to the retarded (ad- 
vanced) Green function. Further we calculate the terms 
of real Pij and imaginary 71 j parts of the integral 
which are defined by the following expressions 



VI. DISORDER EFFECTS. 



/3a (6) 



2 
C2Z? 



3i.T.a(s)ds , 7q(6) 



2 
C2D 



00 

/ 5EQ(s)ds 



A. Multiple core. Clean limit with incoherent 
tunneling. 

We study the disorder effects by introducing the mo- 
mentum scattering first into the tunneling process as de- 
scribed by the incoherent tunneling model. Since the 
tunneling is considered as a perturbation one can assume 
a specular QP scattering at the interface on the bulk side 



and play the role of energy shifting and spectral branch 
broadening, respectively: 



iV,(s,6) 



F7(&)e 



-\s\li2D 



[£- £2(6) -/?(6)]2+ 7^(6) ' 



(85) 



Since parameters /3,7 ~ F/A and t2{b)/T ^ 1 are small 
for b <^ ^2D and |£| > F, the LDOS reaches its bulk value 
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FIG. 8: (Color online) LDOS in logarithmic scale for incoher- 
ent tunneling in the clean limit. Curves, taken for different 
distances p from the vortex center, are vertically shifted for 
clarity. The peaks in LDOS exist up to distances ~ ^2d ■ Here 
A/r = 5, g„ = 1. 



in this region: 



N{p,e) 



Ve2 - r2 



(86) 



Skipping the standard calculations of the above- 
defined integrals ([M)) . we get the final expressions for 
parameters (see Appendix \X\ for details): 



/3 = 






sign(e 4- fib) 



n^r^'n 



Ar 



\Qn 



nb^ 



(87) 



The angular brackets denote averaging over the momen- 
tum Qz along the vortex axis in bulk, fl = deo/db. 
The DOS has a peak of height r/7 at an energy 
e = C2{b) + I3{b) shifted from the standard bound state 
level. This shift resuhs in splitting of the ZBA^^ (Fig. [8]). 
For calculations we use the numerical procedure similar 
to that used earlier for the coherent limit; the induced 
potentials were averaged over the cylindrical Fermi sur- 
face in the bulk. 



B. Multiple core. Dirty SC with clean 2D layer. 

Smearing of the energy dependence of the induced 
potentials caused by disorder becomes even stronger if 
the bulk SC has a short mean free path: £ -^ £,s- In 
dirty limit, the momentum averaged retarded (advanced) 
Green functions are parameterized as follows: 



5f(^)(p)=f3sine^(-^) 



T2 COS 



QBiA)^ 



We put 6^("^) = ±01 -I- 102. The boundary conditions 
for p ^ are 5^^^) -^ ±1, f^^^) j'^^W -^ which 
requires Qi — > 7r/2, 82 — > 0. At large distances 9i — > 
0, tanh82 -^ — e/Aoo for e < Aoo while 9i -^ 7r/2, 
tanh82 -^ -Aoo/e for e > Aoo. Then, 82 = for 
e ^ Aoo, and the Usadel equation becomes^ 



D. 



V28i 



sin(28i) 
2p2 



-2Aosin8i = 



(89) 



FIG. 9: (Color online) The local DOS in logarithmic scale for 
the dirty limit with the parameters A/F = 5, V2f /Vf = 1. 
Curves, taken for different distances p from the vortex center, 
are vertically shifted for clarity. 



The solution of Eq. (|89l) has been found in Ref. ,33> 81 (p) 
monotonously decays from 7r/2 at the origin down to zero 
at yO ^ ^5. The Green functions ([55]) determine the in- 
duced vortex potentials Et — i^Qs- 

For small impact parameter values 6 ^ ^5 we get 
Y}j'^ = and the matching condition takes the form: 



^2D0iso) = 2<(so) / smeds + 2ig{so)b\n[sQ/^s] ■ (90) 

The coefficient C in this case has the only broadened pole 
at e = e2(&): 

C[e-e2(6)+*7]=r, (91) 

where the broadening 



7 



2rx/r2 



hv2F 



sin 8 ds 



where the integral is taken along the trajectory. For |e| < 
r and p < Cs the angle-resolved DOS can be written in 
the form 

r2 7(6)e-^l''l 



N,{s,b) = 



VT^^l^[e-e2{b)]^+j^{b) 



(92) 



Consequently, the LDOS has a peak of the height ~ 
r/7(p) at energy e = 62 (p). 

For the energies above the induced gap e > F and small 
impact parameter values €2{b),'-f{b) <^ T the local DOS 
can be replaced by its bulk value: 



Nip,e) 



V^ 



Y2 



(93) 



For 6 3> ^5 the imaginary part of energy decays expo- 
nentially, and Eq. (|?T|) transforms into Eq. pS)) . 

The numerical results shown in Fig. |9] clearly demon- 
strate the broad peak in the LDOS; this peak shifts and 
becomes sharper as the distance from the vortex center 
increases. For p ^ (,s, the LDOS approaches that ob- 
tained for the clean limit in Figs[6]and|8l For calculations 
we used the standard relaxation method^ of solving the 
Usadel equation in the bulk and the Riccati parametriza- 
tion for Eilenberger equations in the 2D layer. 
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C. Vortex core expansion. Dirty SC and 2D layer. 



VII. DISCUSSION. 



To complete our analysis we briefly discuss the case of 
strong disorder both in the bulk SC and in the 2D layer. 
In this limit our model reduces to the one studied numer- 
ically in Ref. [S^. The condition ^s ^ ^20 — \/^J^2dJ^ 
ensures that the short-distance inhomogeneity in the in- 
duced vortex potentials inside the primary core region 
does not disturb the adiabatic solution based on Eq. (1^ . 
Indeed, for momentum-orientation-averaged Green func- 
tions in 2D layer 



9{P) = 



.92 



^n^ 



-i(p 



he'' 

92 






3(k,r) 



one can derive the equation: 



-2(e + Ei)/2 + 2E2.g2 = 0, (94) 

with E2 = £26""^. This equation is similar to that de- 
rived by Kupriyano"\«^ for a contact of two dirty super- 
conductors. 

Using a standard parametrization g{p) ~ T3 sin 4* + 
T2 cos^e"*"^^"^ and the expressions for the vortex poten- 
tials one can obtain the following equation 



iD' 



2D 



V^* 



sin 2* 



2p2 



2rsin(* - e) - 2iecos* = 



(95) 

where V"^ — p ^dp{pdp) and D2D = fiV2pT/2 - 2D diffu- 
sion coefhcient. Integrating Eq. (|95|) . multiplied by p, in 
a small region around the origin (from p = to the value 
Cs ^ Po ^ ^2d) we find the matching condition for the 
adiabatic Green function (|4T1 142]) : 



D 



2D 



9 T 

dp 



Pa 



Pa 



sin 2^1/ 
2p 



dp 



rPo 
2 / pdp [r sin(* - 9) + ie cos 'J'] 
Jo 



. (96) 



Considering the expansion ^(po) = *o ~ Kp^) with 
K = d'ii{po)/dp ^ ^2D ^^'i assuming vJ/q ^ tt/2 one 
obtains cos *o ~ Pq/{£.2D 1^ (Po/Cs)) ^ 1- This estimate 
confirms the conclusion that the LDOS in the dirty limit 
follows the bulk LDOS pattern scaled with the 2D coher- 
ence length ^2r> to within the second order terms in the 
small parameter po/^2D- 

The resulting problem at low energies e <C Aqo coin- 
cides with that describing a standard vortex in a dirty 
SO^ with the gap value F. Thus, the full disordered 
system should reveal the same LDOS patterns as in the 
bulk case, though scaled with the much larger coherence 
length ^2D instead of ^5 . This vortex-core expansion can 
account for anomalously large vortex images observed in 
MgB2^^ and in high- Tc cuprates^H.. 



The results described above imply that the electronic 
states in the induced superconducting configurations 
strongly depend on the tunneling mechanism and on the 
crystal structure of bulk and 2D materials. The structure 
and symmetry of electronic states can be essentially dif- 
ferent from those in the bulk SC. This imposes severe re- 
strictions on possible realization of various exotic proxim- 
ity electronic states^^"^^ including Majorana states^ and, 
in particular, Majorana states in the vortex cores. Our 
results directly show that the existence of zero-energy 
states in the proximity induced vortex core crucially de- 
pends on the tunneling mechanism underlying the prox- 
imity coupling between the 2D layer and bulk SC. One 
expects that the zero energy core state can exist for co- 
herent tunneling between SC and 2D layer both having 
isotropic Fermi surfaces, provided the symmetry of the 
induced superconducting order permits. 

It is known that a zero energy core state exists for a 
vortex with an odd vorticity in graphene monolayer with 
intrinsic superconductivity^"—. The graphene mono- 
layer with proximity-induced superconductivity thus 
would seems to be a good candidate to look for a zero 
energy state. However, the Fermi surface of graphene is 
highly anisotropic; it lies near the Dirac corners of the 
Brillouin zone with the group velocity directed radially 
from the Dirac points. This group velocity direction does 
not coincide with the direction of the Fermi momentum 
and of the Fermi velocity in the bulk SC as shown in 
Fig. [71 Though the results of the previous sections were 
obtained within the quasiclassical approximation, they 
still can shed a light on the possibility of the zero en- 
ergy state in graphene, especially for sufficient doping 
level when the quasiclassical approximation for graphene 
is justified^'^. In this case the results of Sec lVBI can be 
applied. They show that each state in the induced vortex 
core with energy e is coupled to an infinite set of levels in 
the primary core. It is the integral I odd which accounts 
for these states. Its real part deals with off-resonance 
states with eigen-energies not equal to e, while the imag- 
inary part comes from the resonance state with the same 
eigen-energy e. According to Sec lVBl the real part of the 
integral lodd disappears for e = and 5 = 0. The fate of 
the imaginary part depends on that is the zero energy in 
resonance with any state in the primary core or not. It is 
known that for an s-wave clean bulk superconductor the 
core levels are discrete with a minigap wq ^ ^ jEp and 
no one lies at zero energy. Therefore, if the levels in the 
bulk are not broadened by disorder of by inelastic scatter- 
ing, the imaginary part of lodd does not appear, and the 
zero-energy state seems to be intact. The discrete nature 
of the core states is, of course, beyond the quasiclassical 
approximation. Therefore, the above consideration gives 
only a hint towards the possibility of zero energy state. 
The detailed analysis is needed which would be based on 
the strict quantum mechanical description. Note that an 
alternative possibility to save the zero energy states in- 
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troducing a cylindrical cavity in the bulk superconductor 
has been considered in Refsj^i^. 

Other important feature of induced superconductivity 
in a LD system is an extremely large coherence length 
^2D- It provides a unique possibility to realize vortex 
configurations with quite unusual parameters. Here we 
discuss briefly some configurations which are of interest. 
The detailed analysis of all these situations requires spe- 
cial considerations which we postpone to future work. 
First of all we note that the results of Sections IIIH fVl and 
the following sections are valid for ^2D ^ niin(r„,AL) 
where r„ is the intervortex distance and \l is the Lon- 
don penetration length in bulk SC. If the vortex lattice 
in the bulk SC is dense enough with the intervortex dis- 
tance ^2D "^fy <Si \l, the induced 2D vortex cores may 
start to overlap. The spectrum £2 will then be mod- 
ified due to intervortex tunneling of QPs (see Ref4^). 
The effect of the intervortex QP tunneling should be 
important provided the splitting of the quantized en- 
ergy levels due to this tunneling exceeds the minigap 
value. The splitting can be estimated as rexp[— r^/^2D] 
while the minigap inside the induced vortex core is of 
the order of r^/Hv2Fk2F- Thus, the ratio determin- 
ing the intervortex tunneling efficiency is the exponent 
with a big prefactor, fi.W2Ffc2Fr~^ exp[— r^/^2D]- It is 
this ratio which controls the interplay between the veloc- 
ity of the trajectory precession and QP tunneling speed. 
The changes in the QP spectrum become essential when 
T^v ^ C2D ln(?iW2Ffc2F/r). The minigap in this case should 
vanish according to the analysis in Ref. 0. 

In some cases the 2D coherence length ^20 can ex- 
ceed the London penetration depth A^; this depends on 
the properties of bulk SC and on the tunneling rate F. 
If ^2D , '"u 3> Xl the superconducting velocity vanishes 
along the trajectories with 6 > A^, thus the spectral 
branch £2(6) saturates already for b ^ Xl- 

Our results for coherent tunneling can be directly gen- 
eralized for clean d-wave bulk SCs with isotropic Fermi 
surfaces. However, the incoherent tunneling destroys the 
superconducting coherence in the 2D layer. As a result, 
the branch €2 disappears, while the QP states for e < A 
have finite lifetimes for distances close to the vortex cores 
in bulk SC. 

Considering possible experimental realizations of the 
induced vortex states one has to remember of the finite 
dimensions L of the 2D layer. The large size of the in- 
duced vortex cores can lead to the situation typical for 
mesoscopic superconducting samples when L is close to 
several ^2d's. The criterion when the vortex spectrum 
transformation caused by the boundary effects in such 
systems becomes important can be found using the re- 
sults of Ref. |40. One only needs to replace the gap, the 
coherence length and the minigap by the appropriate val- 
ues in the 2D layer. The criterion appears to be very sim- 
ilar to that describing the efficiency of intervortex tun- 
neling: the mesoscopic fluctuations of quantum levels in 
the 2D core become comparable with the minigap for 

L < 6-D \n{hv2Fk2F/^)- 



In conclusion, the model of proximity coupled 2D layer 
gives the possibility to study theoretically many spa- 
tially inhomogeneous situations including various config- 
urations of induced vortices. Based on this model we 
have presented here description of the vortex core states 
for some typical tunneling mechanisms. In particular, 
our results can be used for interpreting the STM data on 
the vortex LDOS in superconductors through the model 
of a thin proximity layer present at the surface of the 
bulk SC. Effect of a thin non-superconducting proximity 
layer can explain various experimentally observed fea- 
tures of the vortex LDOS and reveals that STM tech- 
nique alone is not sufficient for identifying multicompo- 
nent or anisotropic energy gap. 
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Appendix A: Calculation of self energies for 
incoherent tunneling 

Assuming small impact parameter values & ^ ^5, i.e., 
we calculate in this Appendix the following integrals from 
the main text: 

00 00 

Paib) == 7 — / ^T,ais)ds , 7q(6) = - — / '^T,ais)ds . 

t,2D J k,2D J 



For this purpose we consider the case of the small impact 
parameter values b <S^ ^s'- 



PM 



/: 



2V% f°° I v\\e-^ 



V2F Jo \ '2Q^p'^ 

1-5R- 



v/e2 - n^p^ 



ds , 



where p^ = b^ + s'^. In this case the first term in the 
above integral is determined by s ~ &: 



Tb 



^ir 



-K ' 



Qilp^ 



ds = Tb 



00 



v\i 



\Qnis^ + b^) 

TTVW \ 



ds 



sign(6)r 



2Qn/^ 
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The second one is determined by very small impact pa- 
rameters and reads: 



ds 



vw 



■K 

2 



ds 



where 6g = e^/fi^ — 6^ > 0. As a result, we find: 



2T6d ' 



fii{h) = sign(6)— ( T^X(^'6' - e') 



7/(&) = 



e v\\e 



-K 



V2F Jo \ ^P^ Q^JWif^^ 



x(r!V'-e')) ds 



= sign(6e) 



r' / v\i ,, n\b\ + \e\ 

V2F \ Q^ y^\n%^-e^\ 



Here we use the following expressions for the standard 
integrals: 



V2F XQ" 



Here x(a;) is the Heaviside theta-function, i.e., xi^) = 1 
for a; > and xi^) = for a; < 0. 

After simplifying the expression for /3(6) = f3j{b) — 
Pi{b) we obtain ([57]) . For b > £^s the quantity (3{b) decays 
as exp(— 2&/^5). 

The expressions for imaginary parts hold for any dis- 
tances p because the delta functions in the integrals select 
only the trajectories that pass at small impact parame- 
ters: 



71 W = 



iL r I "II 



-K 



V2F Jo \ Q^n^p^ - £2 




V2F \Q^ ^\n^b^ - £2 



ds , A 

= In ■ 



bo V^^±M ^\n%^-e^\ 



where s^^ax ~ £,s, and 



bo 17W±^{JTb^} " N "" ^\n^b^-e^\ 



The imaginary terms also decay exponentially for b > S^s- 
The expression for 7(6) = 71(6) — 7/(6) gives (|87)) . 
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